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It is shown by numerical solution of Maxwell's equations that, for 
a given wavelength, the degree of confinement of the electromagnetic 
field to the core of a graded-index, single-mode, optical-fiber can be 
optimized by the proper choice of the radial variation of the index. 
Such confinement of the energy to the core helps alleviate loss. The 
fibers considered have zero total dispersion bandwidths in excess of 
100 GHz -Km, at wavelengths between 1.3 pm and 1.55 fim. 

I. INTRODUCTION 

Our earlier work described a method of designing single-mode light- 
guides with zero total dispersion by varying the index profile in the 
core. 1 In the range of wavelengths between 1.3 jum and 1.55 /mi, 
bandwidths in excess of 100 GHz -Km are attainable by balancing 
material dispersion with waveguide dispersion. 2,3 However, one of the 
serious difficulties with single-mode fibers is microbending loss. We 
know that the microbending loss for the case of the step-index, single- 
mode fiber is proportional to A 2 / A 3 . 4,5 ' 6 Here, X refers to the operating 
wavelength and A to the relative index difference which is defined as: 

A = (N COTe - Nclad)/iVc„re . (1) 

Now the design of a single-mode fiber must be such that it prevents 
the field of the fundamental mode (HEn) from extending well into the 
cladding. In other words, the electromagnetic field must be tightly 
confined to the core. To this end, two methods can be considered: one 
is to increase A, which is a common and prevailing method, and the 
other is to change the index profile in the core. This work will focus on 
the latter case by assuming an a index profile, where N = Ni{l — Ar Q ). 
Reducing microbending by profile design might be advantageous be- 

1727 



cause the te and tm modes can be maintained well beyond the cutoff 
point of a step-index, single-mode fiber and the manufacturing toler- 
ances are relaxed. Note that Ni = N COT e and N2 = 2V c iad- 

For a single-mode lightguide having a radially inhomogeneous core, 
it is usually not possible to obtain analytical solutions of closed form 
for Maxwell's equations. Hence, to attain vector electromagnetic field 
distributions of the HEn mode and to determine accurate propagation 
characteristics of a single-mode fiber, we used a numerical method to 
solve the governing equations. 7,8 

II. THEORY 

Our method of solving Maxwell's equations for lightguides has been 
described in our earlier publications. However, we did not consider the 
cladding fields in much detail. For the work to be described in this 
paper, this is essential. Thus, we develop the necessary mathematical 
expressions. 

In an optical fiber having a permittivity e and permeability ju, we 
assume that the outer diameter D is much larger than the core 
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Fig. 1 — A cross-section of a fiber and its cylindrical coordinate system. 
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diameter 2a. In a cylindrical coordinate system, for a position vector 
r having as its components {R, <j>, z), the corresponding components 
of electric and magnetic fields can be written as E = {Er, E& E z ) and 
H = {Hr, H& Hz) (see Fig. 1). However, in obtaining a complete set of 
vector solutions, it is only necessary to find the tangential components 
{E+, E z ) and {H^,, H 2 ) since the radial components Er and Hr are 
linear combinations of the other components. In particular, 



E R = - 
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(2) 
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In addition, the tangential field components are continuous through 
the core-cladding interface and this simplifies the mathematics of the 
boundary value problem. 

In the above equation, Z Q is the wave impedance defined by (/*/e) I/2 
and N is the index of refraction and a function of R. The effective 
refractive index N e is defined by two quantities, /? and k, where /? is 
the propagation constant along the fiber axis and k = 2tt/\\ p is a 
dimensionless quantity defined by Rk. The fundamental HEn mode 
propagates in the fiber when the angular mode number M equals 1. 
Moreover, V < V, must be satisfied. Here Vis, the normalized frequency 
and V c is the cutoff frequency for a single-mode propagation. The V 
value is defined by 



v =*™»JnT=nI. 



(4) 



The input data, a„,„, is the optimum core radius that will give zero 
total dispersion for a given a, A, and A. 

In the most general case, there are two possible solutions to Max- 
well's equations for a guided mode in a lightguide. A general solution 
will be a sum of these two vector solutions. We introduce variable I\ 
to establish the following relationship with the tangential field vectors 
{E+E,)and{Hi,H,}. 
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From eq. (5) we denote the two solutions in cores r,i and r,2. Since 
our earlier work gave a detailed description of the computation of Va 
and r,2, we will avoid repetition of the procedure for the solution in 
the core region. 1 In the cladding region, the two solutions designated 
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by r, 3 and T i4 are given by the following expressions: 7 

N 2 e -K(c) 

T i3 =W 1 U) 



k(c)7 1 (() 

N e 





(6) 



and 



rn-mu). 





iV e 
Y,(l) 

N 2 e -K(c) 



(7) 



(8) 



where k(c) is the dielectric constant in the cladding, and 

£=[2Vl- K (c)] 1/2 p, 
W l U)=K i (Z)/[N 2 e - K (c)] 

rjt) - t-KiW/KiiS) 

where #1 is a modified Bessel function of the second kind and its prime 
denotes differentiation with respect to £. (Note that p, is the value of 
p at the interface.) 

The total solution T can be written in the core and cladding region 
separately. 
In the core region, T is expressed by 

T = AiTa + A 2 Ti2 , (9) 



and in the cladding, T is expressed by 

T = A 3 Ti3 + A 4 T i4 



(10) 



where Aj is an arbitrary constant, j = 1, 2, 3, 4. 

To calculate the field function T, we require basically four input 
data, namely, A, the optimum core radius a op t, N e , and N. Among those 
parameters, calculation of N e has been described in detail in Ref. 1. 
The material dispersion effect is incorporated with Maxwell's equa- 
tions to achieve a high degree of accuracy for N e . This is needed to 
acquire the precise eigenfunctions from eqs. 9 and 10. 

In the design of a single-mode fiber, A is usually specified as an input 
data. It is rather small, ranging from 0.002 to 0.008, since the cladding 
of the fiber is generally made of a high-silica glass. The dispersive 
character of the cladding is well known. 9 Therefore, for a given A, the 
index of the core center N\ can be expressed in terms of N2 by 

*-ciStf 

The dispersive properties of the N 2 in eq. (11) can be described by a 
modified Sellmeier formula. 1,7 
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N 2 = Co + GA 2 + C 2 A 4 + 



+ 



a 



(X 1 - /) (A* - If (A* - /) 



, (12) 



where I = 0.035. The coefficients C, are given in our previous work. 7 

For the index profile, we use a well-known formula that is particu- 
larly useful in fiber design. 



N=Nx 



a 



(13) 



Finally, the dispersion of the index N will be determined by substitut- 
ing Ni in eq. (13) with eqs. (11) and (12). 

At the core-cladding interface, T must satisfy the continuity condi- 
tion of the tangential field components. Consequently, this yields a set 
of simultaneous equations 



AiVu + A 2 T i2 = A 3 T 13 + A4T14 
AiT 2 i + A 2 T 22 = A 3 T 23 + A4T24 
AiT 31 + A 2 T 32 = A 3 r 33 + A4T34 
AiTu + A 2 T i2 = A 3 T 43 + A4T44 



(14) 



To compare the field distributions, it is convenient to introduce the 
following normalized variables into eq. (5). 



E z = 



E z 



E* = 



Er — 



H z = 



Ha — 



H R = 



EM 

_E ± _ 

E40) 

E R 

E H (0) 
H z 
H^O) 

Ha 

HM 

Hr 



(15) 



Hr(0) 



(16) 



III. ELECTROMAGNETIC FIELDS FOR THE HEn MODE IN DISPERSION- 
LESS SINGLE-MODE LIGHTGUIDES 

We begin our study by considering a germania- doped silica light- 
guide with A = 0.002 and A = 1.33 jum. The profile parameters examined 
are a = 100, 2, and 1. Thus, we span the range from rectangular 
through parabolic to triangular index functions. In Table I we calculate 
the values for the radii to make these lightguides dispersion-free. 
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Table I — Radii values for 
dispersion-free lightguides 



100 

2 

1 



4.142 
5.725 
6.294 



The normalized electromagnetic fields as a function of normalized 
radii for these three cases are shown in Figs. 2a, 2b, and 2c. Note that 
in all cases the z fields are much smaller than the other field compo- 
nents. In fact, these fields are less than 3 percent of the tangential 
fields. The magnitude of the R and <£ components of the electric and 
magnetic fields for a given a are all essentially the same. This is not 
likely if the index profile becomes more complex. For example, a profile 
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Fig. 2a— Normalized field distributions of the HE n mode in a single-mode lightguide, 
where A = 0.002, A = 1.33 fim, a = 100, and a opI = 4.142 fim. 
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containing a central "burn out" and ripples, which would be charac- 
teristic of modified chemical vapor deposition (mcvd) profiles, would 
not have such a simple relationship between field components. 

An interesting observation from Fig. 2 is that the slopes of the field 
components at the core-cladding interface change with a. For a = 100, 
the field distribution near the interface forms a cusp, but it rounds 
progressively as a decreases. This is due solely to the index distribution 
in the core of the single-mode fiber. 

Figure 3 shows the normalized transverse components of the elec- 
tromagnetic field as a function of radial distance for the three a values. 
The curves are essentially identical up to R = 3 fim. However, beyond 
that distance they deviate appreciably. Also shown, by the vertical 
lines, are the optimum radii. 
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R/a. NORMALIZED RADIAL DIRECTION 



Fig. 2b — Normalized field distributions of the HEu mode in a single-mode lightguide, 
where A = 0.002, a = 1.33 /im, a = 2, and a op , = 5.725 /im. 
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0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.2 

R/a, NORMALIZED RADIAL DIRECTION 

Fig. 2c — Normalized field distributions of the HE n mode in a single-mode lightguide, 
where A = 0.002, A = 1.33 /mi, a = 1, and a op , = 6.294 fim. 



IV. ENERGY FLOW FOR THE HEn MODE IN DISPERSION LESS SINGLE- 
MODE LIGHTGUIDES 

So far we have only considered the field in the fiber. However, in 
experimental practice it is more convenient to know the field intensity, 
which is the amount of energy flowing through the cross-section of the 
fiber. This can be calculated from the Poynting vector S in Ws/cm 2 . 
The Poynting vector in the z direction, S z is given by, 

S 2 = V 2 (E r H; - E^H%), (17) 



where * indicates the complex conjugate of the variable. 
We define the normalized Poynting vector J by 



(18) 



I=S z /S z (0). 
Figure 4 shows the curves of J versus the normalized radial coordinate 

1734 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1981 



R/a for three different a values. The normalized Poynting vector (field 
intensity) falls off more rapidly with normalized radius for lower a 
values. 10 As in Section III, when a = 100 the Poynting vector develops 
a cusp at the core-cladding interface. We also note a near identity of 
the a = 1 and a = 2 curves. Thus, these two have nearly the same 
focusing power. 

V. DEGREE OF FIELD CONFINEMENT 

We know that the degree of field confinement in a fiber is related to 
its microbending loss. 4,5 In the fabrication of single-mode fiber cables 
for undersea applications, microbending loss has been one of the 
factors that determines the performance of the cable. 
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Fig. 3 — Normalized field components versus radial distance for three different values 
of a. 
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Fig. 4 — Normalized Poynting vector versus normalized radial coordinate for three 
different values of a. 



Figures 3 and 4 indicate that the field or power distribution is largely 
dependent on the index profile. To quantify the focusing power or 
confinement of a lightguide, we introduce a parameter O defined by: 



4> = 



S z RdR 



I 



(19) 



S z RdR 



The parameter <J> represents the degree of power confined to the core 
with respect to the total propagating power. This ratio is plotted in 
Fig. 5, along with a opt , as a function of a. We see that a op t increases 
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with decreasing a and that $ reaches its peak very near a = 2. A 
slightly larger value for <t> occurs if the profile is Gaussian; that is, 



MS)}- 



N ^MMt^IM- 



(20) 



This value of <P is the open circle in Fig. 5. We suspect that this slightly 
larger value may be caused by the close matching of the field with the 
index profile. The Gaussian index profile and the a = 2 profile yield 
~40 percent increase in 4> over the step-index profile case. This may 
help in eliminating microbending loss in single-mode fibers without 
increasing A. 

VI. CORE-TO-CLAD RATIO 

In the design and fabrication of single-mode lightguides, it is custom- 
ary to fix the core- to-clad ratio at 0.1. This value seems appropriate 
for step-index fibers. It is, therefore, quite important to investigate the 
behavior of the evanescent field for graded-index fibers. From eqs. 6 
and 7, we can readily calculate the field intensity in the cladding for 
different values of a and any radius. 

For the cases of a = 100, 2, and 1, including a Gaussian index profile, 
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Fig. 5 — Degree of field confinement in a single-mode fiber versus profile exponent a, 
(solid line). The dotted line shows the optimum core radius corresponding to the a value. 
The open circle indicates the maximum value of <t> obtained from a Gaussian-like index 
profile. 
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Fig. 6— Normalized evanescent field intensity at R » a versus normalized radius for 
three different values of a and a Gaussian index profile. The horizontal dotted line is the 
cutoff level at 7 = 10 7 . 



the results are given in Fig. 6. To compare the core-to-clad ratios for 
different as we define an intensity level / at R :» a equal to IO -7 as 
the cutoff point. This corresponds to R/a ~ 9.3 for a = 100. Accord- 
ingly, Fig. 6 shows that there are substantial differences in those values 
among the four cases. The clad-to-core ratio is reduced to 7.6 from 9.3 
as a decreases to 2. The value for the Gaussian profile is very close to 
that for a = 2. Finally, it is interesting to note that the value of /at the 
core-cladding interface is -0.13 for a = 1 and 2, but it is -0.28 for a 
= 100 (see Fig. 4). 
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VII. VALIDITY OF GAUSSIAN APPROXIMATION FOR THE FIELD 
FUNCTIONS OF HE,, MODE 

As mentioned earlier, it is usually not possible to find an analytical 
expression for the electromagnetic field functions for the HEn mode 
in a single-mode fiber. One exception to this is the step-index profile. 
Therefore, an approximate expression for the field distribution of the 
fundamental mode is frequently used to determine the propagation 
characteristics. A prevailing approximation is a Gaussian-like 
field. 111213 Thus, we can write 



E e = E Q exp 



-61* 

a 



(21) 



GAUSSIAN FIELD 
APPROXIMATION 




--5 



Fig. 7a — Comparison of Gaussian field distribution with exact field solutions. The 
solid lines are the exact values, and the dotted lines are the Gaussian approximation for 
m = 2 and a = 100. 
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where E Q and b are constants. Taking a double logarithm of both sides 
of eq. 21, we can rewrite it as 



where 



lnln(i-) =ln6 + mln(-), 



Em — E M /E . 



(22) 



We introduce the dimensionless quantity U to represent any one of 
transverse electric or magnetic field components, for example, Er. 

Equation (22) was plotted with m = 2 (precisely Gaussian) and 
compared with the exact values. The results are shown in Figs. 7a, b, 
and c for cases a = 100, 2, and 1. The solid lines are the exact values, 
while the dotted lines are from the Gaussian approximation. In all 
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Fig. 7b— Comparison of Gaussian field distribution with exact field solutions. The 
solid lines are the exact values and the dotted lines are the Gaussian approximation for 
m = 2 and a = 2. 
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GAUSSIAN FIELD 
APPROXIMATION 




Fig. 7c — Comparison of Gaussian field distribution with exact field solutions. The 
solid lines are the exact values, and the dotted lines are the Gaussian approximation for 
m = 2 and a = 1. 



cases, the fields of the core region appear in the third quadrant of the 
figures. Those in the cladding region are shown in the first quadrant. 
For the core region, when a = 100 and it is near the core-cladding 
interface, there is satisfactory agreement between the Gaussian field 
function and the exact field function. There is poor agreement near 
the center of the core; this is also evident from Refs. 11 and 12. When 
a = 2 or 1 there is much worse agreement between the exact and 
Gaussian functions in the core regions. For all values of a, the agree- 
ment in the cladding region is extremely poor. It is interesting to note 
that the slope of the exact field functions in the cladding for all a 
values is close to 1. This indicates that the field is decaying exponen- 
tially. 
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VIII. SUMMARY AND CONCLUSIONS 

From a numerical solution to Maxwell's equations, we can accurately 
describe the field distribution of the HEu mode in a single-mode 
lightguide. According to our calculated results, the fraction of power in 
the core reaches its maximum value near a = 2. Thus, for a parabolic 
profile or Gaussian profile, the fraction of power within the core is ~40 
percent larger than that of a step-index core. On the other hand, it is 
clear that the optimum core size increases with decreasing a value. A 
linear index profile (a = 1) provides an optimum core size that is over 
50 percent larger than that of a step-index core. Therefore, in designing 
a single-mode fiber, it is important to remember that one should 
choose a value of a to optimize certain characteristics, such as zero 
total dispersion, te and tm cutoff, core size, manufacturing tolerances, 
field confinement, or microbending loss. 

The work of Marcuse shows that for the case of step-index, single- 
mode fibers, microbending losses can either increase or decrease as a 
function of fiber radius, depending upon the statistics of the axis 
deformation function. 141516 He also concludes that single-mode fibers 
with parabolic-index profiles may have smaller microbending losses 
than single-mode, step-index fibers. If the distortion power spectrum 
peaks sharply at low spatial frequencies, the advantage will be slight, 
if it exists at all. However, for distortions with a wider Fourier spec- 
trum, the parabolic-index fiber should clearly be advantageous. The 
reason for this is that in the case of the distorted parabolic-index fiber, 
the sources of the radiation field are distributed throughout its volume, 
while in the case of the step-index fiber, they are located at the 
waveguide boundary. The constructive interference among the volume 
sources is never as pronounced as among the boundary sources. 

It should be noted that when we increase the mode confinement we 
reduce the field at the core-cladding interface. This reduces the 
strength of the radiation sources because of microbending at this 
boundary. Furthermore, if there is a barrier layer such as B2O3, then 
bulk loss is reduced as well. 

An additional discussion seems in order concerning the important 
te and tm cutoff. As we have previously shown 1 for A = .002 and \ = 
1.33 /tm, the te and tm cutoff is 1.0 /im when a = 100, and it shifts to 
0.85 /im for a = 1. Attempts to increase the field confinement by 
increasing A, while keeping the core radius fixed, may move the cutoff 
rather close to the operating wavelength. Instead of increasing A it 
may sometimes be better to reduce a. The substantial increase in core 
size when a = 2 is an advantage as far as coupling to a source is 
concerned. It has the important added advantage in that the clad-to- 
core ratio is reduced. 

Finally, we must conclude that the Gaussian field approximation is 
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likely to be of value only in the core region of a step-index fiber. It is 
always very poor in the cladding. However, the Gaussian-like approx- 
imation with m 5* 2 may be useful for graded-index, single-mode 
lightguides. For example, for a = 1 and 2 the m values for best fit are 
~1.6 and —1.8, respectively. 
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